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Abstract 

We determine the critical exponents for the XY universahty class in three 
dimensions, which is expected to describe the A-transition in ^He. They 
are obtained from the analysis of high-temperature series computed for a 
two-component A(/)^ model. The parameter A is fixed such that the leading 
corrections to scaling vanish. We obtain = 0.67166(55), 7 = 1.3179(11), 
a = —0.0150(17). These estimates improve previous theoretical determina- 
tions and agree with the more precise experimental results for liquid Helium. 

Keywords: Critical Phenomena, A-transition, High- Temperature Expansion, 
Critical Exponents. 
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The renormalization-group approach to cooperative transitions is one of the most signif- 
icant successes in theoretical physics. It has explained a wide range of phenomena in many 
different fields, ranging from statistical physics to elementary particle physics. It should be 
noted however that the quantitative experimental support for the theory rests primarily on 
its predictions for critical phase transitions. To test the theory one needs the most accurate 
possible values for universal quantities. Here the ^He phase transition is of great utility 
because of the weakness of the singularity of the compressibility of the fluid and of the avail- 
ability of extremely pure samples. Moreover the possibility of performing the experiments 
in space, and therefore in a microgravity environment, reduces the gravity-induced broad- 
ening of the transition. Recently a Space Shuttle experiment performed a very precise 
measurement of the heat capacity of liquid helium to within 2 nK from the A-transition 
obtaining an extremely accurate estimate of the exponent a: 

a = -0.01285 ±0.00038. (1) 

This estimate is extremely precise and represents a challenge for theorists, who, until now, 
have not been able to compute critical exponents at this level of accuracy. 

One of the oldest approaches to the study of critical phenomena is based on high- 
temperature (HT) expansions. In this approach the main hindrance to a precise deter- 
mination of universal quantities is the presence of confluent corrections with non-integer 
exponents. For instance the specific heat is supposed to behave as 

Cp = At"" (l + + Ct . . .) + D + + . . . (2) 

for t = |T — Tc|/Tc — ^ 0, with A 0.5. The presence of non-analytic terms introduces a large 
(and dangerously undetectable) systematic error in the results of the HT series analysis. In 
order to obtain precise estimates of the critical parameters, the approximants of the HT 
series should properly allow for the confluent non-analytic corrections. Integral (also called 
differential) approximants are, in principle, able to describe a behaviour of the type (^ 
(see e.g. Ref. |Q] for a review). However, the extensive numerical work that has been done for 
the Ising model shows that in practice, with the series of moderate length that are available 
today, no analysis is able to predict and take into account non-analytic correction-to-scaling 
terms In order to effectively keep into account these confluent corrections, one 

should use biased approximants, fixing the value of /3c and of the first non-analytic exponent 
A (see e.g. Refs. |U|,0,|TH1). 

To overcome these difficulties, in the early 80s, Chen, Fisher, and Nickel realized the 
importance of studying families of models (specified by some auxiliary parameter) which 
are candidates for belonging to the same universality class. The hope was the possibility 
of locating a parameter value at which the leading non-analytic corrections vanish. If the 
leading non-analytic terms are no longer present, one expects a faster convergence, and 
therefore more precise and reliable estimates of the critical quantities. The method was 
applied to the double- Gaussian and to the Klauder models, both belonging to the Ising 
universality class and depending continuously on a real parameter; it was shown that a 
Hamiltonian for which the leading corrections are suppressed — we will name it "improved" 
Hamiltonian — could indeed be found P,|15],|16|. The crux of the method is the precise 



determination of the optimal value of the parameter appearing in the Hamiltonian. In 
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Refs. [^|T5| the partial differential approximant technique was used; however the errors on 
the improved Hamiltonian were relatively large and the final results represented only a 
modest improvement with respect to standard (and much simpler) analyses using biased 
approximant s. 

In the last few years it has been understood that improved Hamiltonians can be deter- 



mined with high accuracy by means of Monte Carlo simulations [p!7|-pO|. Using finite-size 
scaling methods that are very sensitive to confluent corrections, it has been possible to 
determine various improved Hamiltonians that belong to the Ising universality class, and 
correspondingly precise estimates of critical quantities have been obtained. It is important 
to notice that the estimates obtained using different improved Hamiltonians agree within 
the quoted error bars, confirming the correctness of the error estimates. Similar methods 
have been applied to the determination of critical exponents for dilute polymers in good 
solvents ||2ll . 



The possibility of determining precisely the "improved" value of the parameter for a 



family of Hamiltonians has recently revived the program of Ref . Q . In Ref . we performed 
an extensive analysis for the Ising universality class. Critical exponents and many other 
universal quantities were determined by analyzing HT series for three different improved 
Hamiltonians. For each improved model, we obtained very accurate estimates |^3[ that were 
in good agreement among each other, confirming the correctness of the quoted errors. This 
analysis showed that, once a precise estimate of the improved parameter is available, the 
HT series analysis gives results of a quality comparable with or better than the best Monte 
Carlo simulations. 



As is well known, see e.g. Ref. the A-transition of liquid Helium is expected to be in 
the XY universality class. In order to check if this is really the case, it is important to have 
precise theoretical estimates that can be compared with the experimental results. Recently, 



Hasenbusch and Torok p5| performed a high-precision simulation of the 0(2) 0^ model, 



obtaining an accurate estimate of the exponent v. p = 0.6723 (3) (8). Using the hyperscaling 
relation a = 2 — 3z/, they obtained 

a = -0.0169 ±0.0033. (3) 

This estimate is larger than the experimental result (p, although the difference is barely 
bigger than the error. It is therefore important to further improve the accuracy of the 
theoretical estimates in order to understand whether this small discrepancy is significant. For 



our purposes Ref. is particularly relevant since it provides a quite precise determination 
of an improved Hamiltonian belonging to the XY universality class. In this paper we study 
this improved model and obtain very precise estimates of the critical exponents from the HT 
series analysis. They considerably improve previous HT estimates, showing the effectiveness 
of the approach. In particular we obtain for a 

a = -0.0150 ±0.0017, (4) 

which is in better agreement with the experimental result, although still slightly larger. 
The small discrepancy between the Monte Carlo and the experimental result is significantly 
reduced, providing support to the fact that the A-transition belongs to the same universality 
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TABLE I. Estimates of the critical exponents. See text for the explanation of the symbols 
in the first column. We indicate with an asterisk (*) the estimates that have obtained using the 
hyperscaling relation 2 — a — 3i^ = or the scaHng relation 7 — (2 — 7/)z^ = 0. bcc and sc refer to 
the body-centered cubic and to the simple cubic lattices respectively. 



class of the XY model. The full set of estimates together with other recent results is reported 
in Table |. 

We consider a simple cubic lattice and the Hamiltonian 

n= -(3 E 4.4+^[02 + A(0^-l)2], (5) 

<x,y> X 

where < x,y > labels a lattice link, and (px is a real two-component vector defined on lattice 
sites. Hasenbusch and Torok have performed an extensive simulation of this model and have 
obtained an estimate of the value of A, A*, at which the leading corrections vanish. They 
quote: 

A* = 2.10 ±0.01 ±0.05, (6) 

where the two errors are respectively the statistical and the systematic one. We have con- 
sidered the susceptibility x ^iid the second-moment correlation length ^ defined by 

X = E(^o-0x.), (7) 

X 

e = ^EN'('^o-4). (8) 

6x X 

Using the linked-cluster expansion technique, we have generated the HT expansion of these 
two quantities to 20th order for arbitrary values of A. From the analysis of x C we 
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can obtain directly estimates of 7 and v. In the analysis we used inhomogeneous integral 
approximants (lA). Second- and third-order lA's turned out to be the most stable, especially 



when we biased the approximants |]T6| requiring the presence of two symmetric singularities 



at /3 = ±/5c- This is a natural requirement, since it can be proved that, on bipartite lattices. 



/3 = —/3c is also a singular point associated to the antiferromagnetic critical behavior |^ 
The results we quote were obtained using this type of approximants. 

In Table |I| we give our estimates of the exponents 7 and v that have been obtained 
analyzing the HT series for the improved Hamiltonian (|^) with A = A*. We quote two 
errors: the first one is related to the spread of the approximants, while the second one 
gives the variation of the exponent when A varies between 2.04 and 2.16, cf. Eq. (|^). It 
should be noted that the first error is somewhat larger than the second one, and therefore 
it is important to further extend the HT series. However, since the second error is far 
from negligible, a substantial reduction of the uncertainty also requires a more accurate 
determination of A*. Note that this might also reduce the first error, since it would enable 
us to work closer to the exactly improved Hamiltonian. We do not present details on the 
generation and analysis of the HT series and we refer the reader to App. A of Ref. |^ . 



From the estimate of v we can obtain the exponent a assuming the validity of the 
hyperscaling relation 

a = 2 - 3z/. (9) 

In principle it should be possible to determine a directly from the specific heat, or from the 
singularity of the susceptibility at the antiferromagnetic point |]3^. We tried this second 
method obtaining only a rough estimate: a = —0.02(2). 

From the estimates of 7 and u it is possible to obtain the exponent rj, using 

7 = K2-^). (10) 

However, it is not clear how to set the error bar on the result. One can use the independent- 
error formula taking into account the error on 7 and u, but this may be an overestimate 
since 7 and u are correlated. To obtain an estimate of rj with a smaller, yet reliable, error 
bar, we used the so-called critical-point renormalization method (see e.g. [@] and references 
therein). The value that is quoted in Table |I| has been obtained with this method. It is 
compatible with the estimate obtained using the scaling relation (p!0D , but it has a smaller 
error bar. 

In Table | we report a summary of the most precise estimates that have been obtained 
in the last years. When only u 01 a was reported, we used the relation (^ to obtain 
the missing exponent. Analogously if only 77 or 7 was quoted, the second exponent was 
obtained using the scaling relation ([lOD ; in this case the uncertainty was obtained using 
the independent-error formula. The results we quote have been obtained from the analysis 
of the HT series of the XY model (HT), by Monte Carlo simulations (MC) or by field- 
theory methods. The field-theory results have been derived by resumming the perturbative 
expansion in fixed dimension d = 3 ((^-expansion) , or the the expansion in e = 4 — d. For the 
e-expansion we quote two numbers, corresponding to an unconstrained analysis ("free"), and 
to a constrained analysis ("be") in which the two-dimensional values of the exponents are 
taken into account. Our final results for 7 and u improve the existing theoretical estimates. 
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They are somewhat lower than previous HT results but they are in full agreement with the 
field-theory and the most recent MC estimates, as well as with the experimental results for 
the A-transition, showing the expected universality. 
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